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Abstract. The nodal domains of eigenvectors of the discrete Schrodinger operator 
on simple, finite and connected graphs are considered. Courant's well known nodal 
domain theorem applies in the present case, and sets an upper bound to the number of 
nodal domains of eigenvectors: Arranging the spectrum as a non decreasing sequence, 
and denoting by i/„ the number of nodal domains of the n'th eigenvector, Courant's 

■ theorem guarantees that the nodal deficiency n — Vn is non negative. (The above 
\ applies for generic eigenvectors. Special care should be exercised for eigenvectors with 

vanishing components.) The main result of the present work is that the nodal defi- 
ciency for generic eigenvectors equals to a Morse index of an energy functional whose 
value at its relevant critical points coincides with the eigenvalue. The association of 
the nodal deficiency to the stability of an energy functional at its critical points was 

■ recently discussed in the context of quantum graphs [T] and Dirichlet Laplacian in 
bounded domains in [S]. The present work adapts this result to the discrete case. 

■ The definition of the energy functional in the discrete case requires a special setting, 
"j^ \ substantially different from the one used in [1] |^ and it is presented here in detail. 



o 
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^ ■ 1. Introduction 

^ \ Courant's nodal domain theorem can be viewed as a generalization of Sturm's os- 

\ dilation theorem to Laplace-Beltrami operators in higher dimensions. In the one- 

dimensional case, after ordering the spectrum of the Sturm-Liouville operator as an 
increasing sequence, the oscillation theorem guarantees that the n'th eigenfunction /^"^ 
flips sign = n — 1 times in the open interval. Equivalently, the number of nodal do- 
mains Un — defined as the number of intervals where /^"^ is of constant sign — equals 
^ ■ n. The Sturm oscillation theorem can be written concisely as 

■ (1) n = l'n = (pn + l ■ 

In higher dimensions, the two equalities in ([1]) have to be modified. For d > 1 there is 
no natural analogue of (/)„, and therefore the right equality has to be discarded. Courant 
[TU] showed that the left equality cannot hold in general, and it should be replaced by 
a bound: 



(2) n > z/, 



where Vn is the number of maximal connected components of the domain on which the 
n'th eigenfunction has constant sign. Later studies [21^ ,22j have shown that equality 
holds only for finitely many eigenvectors. These states are referred to as Courant sharp. 

Courant's nodal domain theorem was extended to the Laplace operator on metric 
(quantum) graphs [13j and on discrete graphs [11] (see also references therein). The 
latter are the subject of the present work, and they will be discussed in detail in the 
next section. 
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The interest in counting the nodal domains increased in the mathematical and phys- 
ical communities when it was realized that the nodal sequence {i'n}'i^=i stores met- 
ric information about the manifold where the Laplace-Beltrami operator is defined 
[3 El El El [ini [IT]. It was shown, in particular, that in some cases the nodal sequences 
of isospectral domains are different [T71 HI [3] and that in some other examples one 
can uniquely reconstruct the domain geometry from the given nodal sequence |18) . 

A new point of view was introduced in the pioneering article of Helffer, Hoffmann- 
Ostenhof and Terracini [T5|, where a variational approach was used to locate the 
Courant sharp eigenfunctions in the spectrum. Helffer et. al. investigated the Dirichlet 
Laplacian in a bounded domain 1] G M'^ with d > 2. Partitioning Q arbitrarily into 
u sub-domains Q^, they studied the lowest Dirichlet eigenvalue Ai(f2fc) for each of the 
sub-domains. The maximal value amongst the Xi{Qk) for a given partition V (denoted 
as A(z/;'P)) can be viewed as the "energy" of the partition. Helffer et. al. proved that 
the partitions that minimize A(z/; V) coincide with the nodal partition induced by a 
Courant sharp eigenfunction if and only if the minimizing partition is bipartite. 

We note that that the restriction of an eigenfunction /"^"^ to a nodal domain is the 
ground state of the Dirichlet Laplacian on fi„. The corresponding ground energy Ai(fifc) 
is equal to the eigenvalue A„(f2) on the entire domain. Thus a partition corresponding 
to an eigenfunction will have a special property: Ai(fifc) is the same for all k. A partition 
with this property will be called an equipartition. 

In two recent papers [T, ^ the approach of Helffer et. al. was broadened in a substan- 
tial way. It was shown that the functional A(z/; V), when restricted to a submanifold of 
equipartitions, becomes smooth. One can then study the critical points of A(z/;P), i.e. 
the points where the variation with respect to perturbations of the partition boundaries 
vanishes. It was shown that the critical z/-partitions that are bipartite are in one-to- 
one correspondence with the nodal partitions of eigenfunctions with u nodal domains. 
Moreover, the Morse index /i„ at the critical partitions equals the nodal deficiency of 
the corresponding eigenstate, 

(3) fin = n-Un . 

Thus, in the space of bipartite equipartitions, the critical points in the landscape of 
A{iy;V) are at eigenfunctions, and their stability (number of directions at which the 
critical point is a maximum) determines the nodal deficiency. The result was proved 
for the Laplacian on metric graphs (quantum graphs) in [Jj and was then shown to hold 
on domains in R'^ in [6j. The present work complements the above mentioned articles 
by showing that ([3|) applies also for the discrete Schrodinger operator on finite graphs. 

The variational approach used in \V, '6] crucially depends on the ability to smoothly 
change the boundaries of the domains. There is no direct analogue of this for discrete 
graphs. In the present work we show how one can use local variation of the potential 
in place of the variation of the partition boundaries. This is the chief new element 
introduced in this article, and it enables us to arrive at the main result of the present 
work, namely equality (jSj) for generic eigenvectors of the graph Shrodinger operator. 
The proof is provided in theorem 14.31 

There is another special feature that distinguishes graphs (metric or discrete) from 
other manifolds: one can describe the nodal structures in terms of points where the wave 
function changes its sign (for discrete graphs these will correspond to changes of the sign 
of the vertex eigenvector across a connecting edge). This allows us to reinterpret the 
bound on the number of nodal domains as a generalization of Sturm theorem. Indeed, 
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we show that the right equahty in ([T]) is replaced by upper and lower bounds on 0„, 
given in equation (129|) , Theorem 13.41 

The paper is organized in the following way. The next section provides a few defini- 
tions and known facts from spectral graph theory which are necessary for the ensuing 
discussion. The construction of partitions and the analogue of boundary variations 
require some ground work which is carried out in the section [3l The main results will 
be formulated and proved in section |H 

2. Definitions and general background 

In the present chapter we provide a few definitions and facts which set the stage for 
the subsequent discussion. 

A graph G = {V,S) consists of a set of vertices V = V(G) and a set of connecting 
edges g = £{G). We shall use: V{G) = \V{G)\ and E{G) = \S{G)\. An edge e G £{G) 
connecting the vertices i,j G V{G) will be also denoted by e = = We 

will use the notation z ~ j to indicate that the vertices i and j are connected (are 
neighbors). A graph is simple if no more than a single edge connects two vertices and 
no vertex is connected to itself; otherwise the graph is a multi-graph. A graph is said 
to be connected if there exists a path of connected vertices between any two vertices 
in V. Unless otherwise specified, the graphs we consider here are finite, connected and 
simple. 

The connectivity of a graph G is summarized by the adjacency matrix A{G). 

l-^fC'l-^' = { I: othe^iL 

The degree of the vertex i is the number of its neighbors; it can be expressed via the 
adjacency matrix as 

The (first) Betty number of a connected graph is defined as 
(5) PiG) = E{G) - ViG) + 1 . 

This is the number of independent cycles on the graph. 

A subgraph G" C G is itself a graph, defined by a subset of the vertex set V(G") C 
V{G) and a subset of the edge set S{G') C £{G). A factor (or spanning graph) of G is 
a subgraph, G' C G, which shares with G the vertex set, V{G') = V{G). By removing 
P{G) edges one can generate a factor which is a tree (a spanning tree). There exist 
more than one spanning tree. In what follows we shall often construct a sequence of 
connected factors by starting with G and removing edges one at a time, while keeping 
the resulting factor connected and ending with a spanning tree after /3(G) steps. 

A class of factors which will play a prominent role here are zz-partitions. 

Definition 2.1. A v-partition of the graph G, denoted by V, is a factor consisting of 
V disjoint subgraphs, P^, such that for any two vertices i and j in the same connected 
component Pjt, if («, j) G £{G) then they are also connected in P^. In other words, 
only the edges running between components were removed from G to form V. We write 

We associate to any partition V a multi-graph G{V) with u vertices. The vertices of 
G{V) are in one-to-one correspondence with the disjoint subgraphs of V, while the edges 
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of G{V) correspond to the edges that were removed from G. A partition is bipartite if 
its graph G{V) is bipartite, or a tree if G{V) is a tree, etc. Note that if G{V) is a tree, 
each pair of nodal domains is connected by at most one edge. Edges of this kind are 
referred to as bridges in the graph theory hterature. 

2.1. The Schrodinger Operator. The discrete Schrodinger operator acts on the 
Hilbert space of real vectors / G R^^'^^ where the components are enumerated by the 
vertex indices i. The Schrodinger operator is a sum of two operators - the Laplacian 
and an on-site potential. The Laplacian is usually defined as 

(6) L{G) = -A{G) + D{G) 

where A{G) is the adjacency matrix of G and D(G) is a diagonal matrix with [Z)(G')]jj = 
di. 

The on-site potential, Q{G) is a diagonal matrix [Q{Gy\ii = G M. However, since 
both D and Q are diagonal matrices and we will not impose any restrictions on the 
choice of site potentials g^, we can "absorb" the degree matrix in the potential and define 
the Schrodinger operator 

(7) H{G) = -A{G) + Q{G). 

This convention significantly simplifies the notation later. The action of the Schrodinger 
operator (which we will also call the Hamiltonian) on a vector / is given by 

(8) [H{G)f\ = -Y^f, + q,h. 

The spectrum of ([7j) is discrete and finite. We denote it as 

aiG) = {A„(G)}ri? with Ai(G') < X^iG) <...< A„(G). 

The eigenvector corresponding to A„(G') will be denoted as /*^") = {fj;^^}]^J:f\ 

Remark 2.2. The definitions above can be generalized by associating positive weights 
rriij = rrij^i to the connected bonds in G. The weighted adjacency matrix is defined as 
[M{G)]ij = [A{G)]ijmij . The results derived in the present paper are valid for this 
generalized version of the Schrodinger operator. However to simplify the notation we 
shall present the results for the case of uniform unit weights. 

2.2. Nodal Domains. Let /:(?—)■ R^^'^), / = {fi}]^^^ be an arbitrary real function 
on the graph G. We define a strong nodal domain as a maximally connected subgraph 
of G such that on all of its vertices the components of / have the same sign. Vertices 
where / vanishes do not belong to a strong nodal domain. A weak positive (respectively 
negative) nodal domain is defined as a maximally connected subgraph Y such that for 
all i G V{Y), fi > (respectively fi < 0). In this weak nodal domain must 
contain at least one vertex where / is non-zero. We will seek to understand the number 
of nodal domains of n-th eigenfunction of the Schrodinger operator ([7]). 

Definition 2.3. We will call an eigenfunction non- degenerate if it corresponds to a 
simple eigenvalue and its values at vertices are non-zero. 

In this paper only non-degenerate eigenfunctions / will be considered^ In such 
cases the strong and weak nodal domains are the same and we denote by //(G; /) their 
number. Thus, any f : G ^ (]R\0)^''*^^ induces a bipartite iy{G; /)-partition of G with 



This behavior is generic with respect to perturbations of the potential Q. 
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components which are the nodal domains. The number of edges of G which are deleted 
to generate the partition will be denoted by /), it is the number of sign flips of 
/ that occur along edges of G. The number of independent cycles in G where / has a 
constant sign (i.e. cycles that are contained in a single nodal domain) will be denoted 
by £(G; /). The following identity relates these quantities [2]: 

(9) v{G- f) = /) - /3(G) + i{G; /) + 1 . 

We shall focus on partitions induced by eigenvectors /("^ of ([7]). Because of the 
special role played by this function, we shall use the abbreviations i/„ = ^{G; /^"■'), (f)n = 
(j){G] Z*-"^) and = £(G; Z*-"^). Here, Courant's theorem [11], and its extensions [5] state 
that the number of nodal domains is bounded by 

(10) n - f3{G) <Un<n . 

Note that the lower bound on the number of nodal domains is not optimal. For well 
connected graphs such as ci-regular graphs, fi{G) = -)- l and hence, for c? > 4, 

we have n — f3{G) < for all n < V{G). An improved lower bound will be derived 
below. 

Following [11 [15], we define. 
Definition 2.4. An eigenvector of the Hamiltonian ([7j) is Courant sharp if z/„ = n. 



The eigenvectors for the discrete Hamiltonian on trees are all Courant sharp [12] (see 
also |5j). While it is a special case of ( fTOj) with /3(G) = 0, it is the first step in the proof 
of the lower bound in ( fTO|) for general /3. 

The chromatic number x of graph can be used to give a bound on the number of 
nodal domains [20j : < F — x + 2. Hence, the only graphs for which the highest 
eigenvector (with n = V{G)) can be Courant sharp, are bipartite graphs where x = 2. 

3. Edge manipulations and the parameter dependent Hamiltonian 

The main result of this paper is stated in Theorem 14. 3[ Its formulation requires 
definitions, concepts and facts which are provided in the first half of this section. To 
gain some intuition to the more formal discussions, we start by making the following 
simple observation. 

Let / = be an non-degenerate eigenvector of the Hamiltonian H{G), see Defi- 
nition [2]3l Let G' be a connected factor of G obtained by deleting an edge e = 
Following [5j, we will modify the potential at the vertices i and j in such a way that 
/ is an eigenfunction of the factor G' with an eigenvalue that equals A„(G). That is, 
there is m, such that Am(G') = A„(G). Note that n and m are not necessarily equal, 
and the other eigenvalues and eigenf unctions of H{G') need not coincide with those of 
H{G). To work out the necessary modification to the potentials we rewrite ([7]) at site i 

(/7(G)/), = -5^/, + g,/, = - ^+U-J)f^ 

k^i kr^i, k^j ^ •^^ ^ 

and similarly at site j. 



[H{G)f), = - Y h+U-T)f^ 



k^j, kj^i 

Thus, 

(11) H{G)f=(^-A{G') + Q{G'))f 
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where the potential Q{G') coincides with the original potential Q{G') on all vertices 
except the vertices i and j where the potential is modified to be 

(12) qi = qi-Y = ii ~ T- 

The operator which multiplies / on the right hand side of ( JTTj) above is a Hamiltonian 
operator for the factor G', 

(13) H{G') = -A{G') + QiG') . 

Clearly / is an eigenvector of H{G') corresponding to an eigenvalue which equals Xn{G) 
but whose position m in cr{G') is not necessarily n. 

This formal exercise acquires more substance once the modified potentials are defined 
by replacing in (IT^ the quotient y- with a real non zero parameter a and y- with 1/a. 

The resulting parameter dependent Hamiltonian H{G'; a) for the factor G' is defined 
as 

(14) H{G';a) = H{G) + B{a). 

where the perturbation B = B{a) has only four non-vanishing entries, 

(15) 5,,, = -a, B,,, = 1 

(16) B,, = 1 B,j = -1/a 

The matrix B{ij) has rank 1, with a single non vanishing eigenvalue —a — ^. 

The following theorem is due to Weyl (see, for example, [T6|) and will be extensively 
used in the sequel. 

Theorem 3.1. Let A,B,G be real and symmetric N x N matrices with G = A + B . 
Denote by ai < a2 < • ■ ■ and ci < C2 < ■ ■ ■ cn the spectra of A and G respectively. 
If B is of rank 1 then the spectra of G and A interlace in the following way: 

(17) ■ ■ ■ < ak-i < Ck-i < ak < Ck < ■ ■ ■ if B >0, 

and 

(18) • ■ ■ < ak-i <Ck<ak< Cfc+i < ■ ■ ■ if B < 0, 

After all these preparations we can arrive to the following theorem. 

Theorem 3.2. Let G be a simple, connected graph with a Hamiltonian H{G), and let G' 
be a connected factor obtained by deleting the edge {i,j) with parametrized Hamiltonian 
H(G';a) defined by (fT^ -(fT6 |) . Let Xm{G';a) and f'^"^\a) be the m'th eigenvalue and 
the corresponding eigenvector of H{G';a). Consider the the eigenvalue Xm{G';a) as a 
function of a. Its critical points a~^ and a^ satisfy the equation 

(19) a^^^' ^ ' 



with the corresponding sign. 

Furthermore, as long as a+ ^ and finite, Am(G'';a^) is an eigenvalue of H{G) 
whose position in a{G) is n G {m — l,m,m + 1}. The eigenvector f^"^\a^) is the 
corresponding eigenvector of H{G). 

Conversely, if f^""^ is an non- degenerate eigenvector of H(G) then 

An) 
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is a critical point of \m{G'\ a) where m is such that n G {m — 1, m, m + 1}. 

Proof. If f^'^\a) is normalized, by first order perturbation theory we have 

(21) <l\m{G',a) ^ / fi^)[a), ^^(G"; «) j(m)^^- 

da \ ' da 

The only a dependent entries of H[G'] a) are the off-diagonal terms in B, see (fT5|) -(|T6 
Hence 

2 



(22) da = - 

which is equal to zero at the solutions of equation ( !T9|) . 

From now on, we consider only the relevant critical points a~^. At one can easily 
check that 

5(a+)/(™)(a+) = 

and therefore 



XUG'; a+)f^"'\a+) = H{G'; a^)f^^\a+) = H{G)f^"'^i 



a 



Conversely, if /^"^ is an non-degenerate eigenvector of H{G) then it was shown in 
equation (fTT]) that /^"'•' is an eigenvector of H{G'] a) for a determined by equation (l20l) . 

Finally, the matrix B is of rank-1 and Weyl's theorem 13.11 guarantees that the spectra 
of H{G) and H{G') interlace according to ( fTT|) or ( fT8l) . showing that n G {m — 1, m, m + 
1}. □ 

Further insight can be gained by examining closely how the position of the eigenvalue 
in the spectrum and the number of nodal domains change depending on the sign and 
the type (minimum or maximum) of the critical point. Note that we will shift the point 
of view and will use G' as the base graph, i.e. we will investigate the change in the 
above quantities when an edge is added back. 

Let a^ denote a relevant critical point of \m{G'\ a). We start with the case a'^ > 0. 
Weyl's theorem (fT7|) implies that for all a > (that is S < 0), 

At the critical point a = a'^ the value Xm{G'; a^) belongs to the spectrum of H{G) and 
therefore equals either Am-i(G) or Xm{G). Hence, Xm{G';a'^) is a maximal value of 
Xm{G'] a) if it equals Xm{G) and a minimal value if it equals Xm-i{G). For future use 
it is convenient to introduce the following notation. Let M{ij) take the value +1 if the 
critical point is a maximum, and if it is a minimum. Let An{ij) stand for the shift 
in the position of the eigenvalue in the spectrum, from its position in cr{G', a'^) to its 
position in o"((j>'). One can summarize the findings so far by the following statement: 

(23) For a'^ > : M{ij) - An{ij) = 1. 
Similarly, for a'^ < 0, Weyl's theorem fllSp implies that 

A™(G)<A^(G";a)<A„+i(G). 

Hence, Xm{G'] a'^) must attain either the maximal value Xm+i{G) or the minimal value 
XmiG). Using the same notation as above we find: 

(24) For a^ < : M{ij) - An{ij) = 

More information on the transition from G to G' is gained by viewing the eigenvector 
/(")(«'=) first as an eigenvector on G' and then as an eigenvector on G. This dual view 
point is now applied to follow the variation in the number of cycles with constant sign 
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and the number of nodal domains, as one counts them with respect to G' or to G. It is 
important to remember that the sign of a'^ is the relative sign of the i and j components 
f^"^\a^) across the edge (i, j). 

Starting with < 0, the number of nodal domains and the number of loops of 
constant sign are not affected by the transition from G' to G since fj/fi < implies 
that i and j are in nodal domains with different signs. Denoting by Ai{ij) the change 
in the number of loops with constant sign, and by Az/(zj) the change in the number of 
nodal domain, we have 

(25) For a'' < : Ai{ij) - Au{ij) = . 

On the other hand, for a'^ > either the number of nodal domains or the number of 
cycles of constant sign will change. Indeed, the edge either connects two vertices 
that already belong to the same nodal domain, increasing 

by 1 (and leaving 

iy{f^"^^) unchanged) or it connects two nodal domains of the same sign, in which case 
z/(/*^™^) decreases by 1 while i(f^™-^) remains constant. This leads to 

(26) For a^>0 : Ai{ij) - Au{ij) = 1 . 

Comparing ([21]) to (125 p and (125]) to ( 12B]) . the observations above can summarized by: 

(27) Ae{tj) - Au{ij) = M{tj) - An{tj) , 

which is valid for both signs of a'^. 

The discussion so far centered on a factor G' which differs from the original graph 
G by the deletion of a single edge. However, by successive applications of Theorem 13.21 
it can be generalized to any connected factor of G obtained by the elimination of an 
arbitrary number of edges while modifying the appropriate vertex potentials. The set 
of parameters will be denoted by a = {oe '■ e G S{G)\S{G')}, the parameter dependent 
Hamiltonian is H{G', a), and Xm{G'; a) is its m'th eigenvalue. 

Corollary 3.3. Let G be a graph as previously and G' a connected factor obtained by 
deleting k < P{G) edges from G. Let a'^ be a relevant critical point of Xm{G'] a). Then, 
provided that none of the components of f^'^\a'^) vanishes at vertices where edges were 
added, \m{G'\a'^) = A„(G) G cr{G) with \n — m\ < k. The corresponding eigenvectors 
are the same. 

The lower bound on the number of nodal domains, the left part of equation (|T0]) . 
was proved by chaining the operations of edge removal |5]. In fact, a more careful 
book-keeping allows us to sharpen the lower bound (see jl] for the same inequality in 
the context of quantum graphs). 

Theorem 3.4. Let /("^ be the n-th eigenfunction of the Hamiltonian H(G) such that 
the corresponding eigenvalue is simple and f^'"''^ has no zero components. Then the 
number of nodal domains of G with respect to f^^^ satisfies 

(28) n - {(3iG) - Q <vn<n. 

Correspondingly, the number 0„ of edges across which the eigenvector f^^^ changes its 
sign satisfies the bound 

(29) n-l<4>n<n + {f3{G) - Q - 1- 

Remark 3.5. Note that the lower bound on the quantity n — {P{G) — £„) is always 
non-negative. 
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Remark 3.6. Equation (p9|) is the generalization of Sturm's oscillation theorem to dis- 
crete graphs. 



Proof of Theorem 3.4 We cut /3 edges of the graph G, modifying the potential ac- 
cordingly, until we arrive to a tree T, such that / = /*^") is it eigenfunction. It is 
eigenfunction number m and, since it is a tree, Fiedler theorem [12j implies it has 
Vrn{T) = m nodal domains. Also, irn{T) = 0, because there are no cycles on a tree. 
We rewrite equation (l27j) in the form 

Az/(e) = An{e) + A£{e) - M(e), 

where e is the removed edge. Adding back the removed edges one by one and adding 
the above identities to the equation //^(T) = m we arrive at 



Since the number of maxima in the sequence is at most /3, the number of nodal domains 
I'niG) is at least n + — /3, proving inequality (125]) (the upper bound is due to [11]). 
Substituting equation for we obtain inequality fl29p . □ 

Corollary 3.7. If f^""^ satisfies the conditions of Theorem \3.4\ and its nodal partition 
graph G is a tree (or, equivalently , the edges on which f^"'^ changes sign do not lie on 
cycles of the graph), then /'•"^ is Courant- sharp: z/„ = n. 

Proof. Indeed, if G is a tree, then no cycles are broken when removing the edges con- 
necting the nodal domains and = □ 



4. Critical Partitions - the main theorems and proofs 

So far we discussed the reduction of a graph to its connected factors. The generation 
of partitions (disconnected factors) requires the introduction of some more concepts 
and definitions. Let P be a bipartite z/(G)-partition of G, see Definition 12.11 and the 
discussion below it. We denote by £{V) the edge set of V. Since we wish to charac- 
terize the partitions that appear as nodal domains, we will ofter refer to the connected 
components Pk of V as "domains". 

Let (f){G;V) denote the number of edges in S{G)\S{V). Construct the Hamiltonian 
H{V] a), a E W^^^^ where for each deleted edge e = i < j, the potentials ae and 

are added to the vertices i and j as was explained above. Note that the ordering of 
the vertices in an edge can be chosen arbitrarily; we chose i < j for definiteness. The 
Hamiltonian H{V] a) is block-diagonal, namely 

H{V;a) = ^H{Pk;a). 

k=l 

Let Xi{Pk;a) be the lowest eigenvalue of H{Pk;a) with the corresponding eigenvector 
denoted by g{Pk', a). Note that since g{Pk', ce) is a ground state, all its components have 
the same sign. We extend g{Pk', a) to all vertices of the graph G by setting gs{Pk', a) = 
for s ^ V(Pfe). We assume g{Pk', a) is normalized and non-negative, 

(30) Yl \gsiPk;a)f = J2 \gs{Pk;a)\' = l, gs{Pk;a)>(}. 

seV(Pfc) seV{G) 
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Definition 4.1. Let V he a bipartite z/-partition of G witli tlie parameter dependent 
Hamiltonian H{V; a). An equipartition is a pair {V, a'), where a' G M'^('^''^) is a vector 
of parameters such that all the lowest eigenvalues Xi{Pk', a') are equal, 

(31) Ai(Pi; a') = X.iP^; a') = ■ ■ ■ = X,{P,; a'). 

Consider the set Q(V) C M'^('^'^) in the space of parameters a where equation (13T]) 
is satisfied (equivalently, all a such that {V,a) is an equipartition). On this set, the 
function A{V; a') = Xi{Pk', «') will be referred to as the equipartition energy. Obviously, 
the index k is arbitrary. 

Intuitively, the {u — 1) equalities in (13T]) reduce the number of independent variables 
required to define the equipartition energy from (f){G;V) to t] := (j){G;V) — {v — 1). 
Notice that, because of ([H]), 

(32) r/ = 0(6-; V) - {v - I) = I3{G) - £{G; V) > . 

If the partition V is induced by an eigenvector f^^^ of H{G) then it also generates 
an equipartition with the vector of parameters defined by 

An) 

(33) "e = ^, for each e = (z,j)G^(G')\^(P), i< j. 

f i 

Then the equipartition energy coincides with the eigenvalue of G, A(P; a'^) = Xn{G). 
In the sequel we shall consider the equipartition energy in the neighborhoods of these 
special points. Note that in the definition above the vertices i and j belong to different 
nodal domains. Therefore the special values of the parameters Ue are negative. This is 
why from now on we restrict our attention to the negative subspace of the parameter 
space, 

The main results of the paper can now be formulated in terms of two theorems which 
distinguish between the cases = and r/ > 0. 

Theorem 4.2. Let V he a bipartite v-partition of G with t] = (f){G\V) — (z/ — 1) = 0. 
Then V is a tree partition and there exists at most one equipartition (P, 7), such that 

7 G M!^'''^'^''. The value X = A(P;7) is the v-th eigenvalue in the spectrum of G, and 
the corresponding eigenvector is Courant sharp. 

Proof. By definition of the partition graph G{V), it has u vertices and (j){G]V) edges. 
Therefore its Betti number is (f){G; V) — u +1 = and it is a tree (see an illustration in 
Fig. m^b) ignoring the dashed edges). 

Let 7 be a vector of negative parameters such that {V, 7) is an equipartition. The 
functions g{Pk,'y), normalized as in ( 130|) . can be used to construct an eigenvector of 
H{G). Consider 

(34) f = J2tkg{Pk;i), 

q=l 

where the coefficients tk are determined as follows. Pick arbitrarily one of the domain 
in V, say Pi, to be the "root" of the partition tree and let ti = 1. For every domain Pr 
adjacent to Pi on the partition tree, let e = {i,j) be the unique edge connecting Pi to 
Pj. with i & Pi and j & Pr. Without loss of generality assume that i < j and take 

, , 9i 

tr = tile — - 
9j 
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Continue the process recursively to domains at larger distances from Pi. The tree 
structure guarantees that this covers the entire set of domains on the graph without 
ever reaching a domain for which the variable t has been previously computed. Note 
that the values of t alternate in sign, leading to the function / whose nodal partition 
is precisely V. We will now show that the resulting function /, equation (15^ . is an 
eigenvector of H{G). 

Start with the Hamiltonian H{V] 7) and observe that since / is composed out of 
eigenvectors of the connected components of H{V] 7) that share the same eigenvalue, 

if(P;7)/ = A(P;7)/. 

We will now add the edges from £{G) \ £{1^) one by one. For simplicity of notation 
we will consider the addition of the edge e = (i, j) between the components Pi and Pr- 
The addition of the edge results in a matrix B defined by (fT5l) - (fT6|) with a = 7e. The 
relevant elements of the vector / are 

fi = hgi, fj = Ugj = ti-feQi- 

A direct computation shows that Bf = and therefore / remains an eigenfunction 
of the modified graph with the same eigenvalues. Proceeding similarly with the other 
edges we conclude that H{G)f = A(V; 7)/. 

The function / has u nodal domains with respect to the graph G. From Corollary 13. 71 
we get that the function / as an eigenfunction of H{G) is Cour ant-sharp, that is, its 
index is v. Finally, since there is at most one Courant-sharp eigenfunction with v 
domains, the equipartition is unique. □ 

Let Z*^"-* be a non-degenerate eigenvector corresponding to the eigenvalue A„(G'). If 
the corresponding partition V is not a tree, there will be other equipartitions locally 
around the special point (133|) . Namely, we shall consider a ball B^{a'^) of radius e 
centered at defined by (133|) . The value of e > is chosen sufficiently small so that 
the following two conditions are satisfied: (z) the variation in the equipartition energy 
A(P; a) is smaller than the minimal separation between successive eigenvalues in (j{G) 
and [a) none of the hyper-planes = intersect the ball. The discussion which follows 
is restricted to a G B^{a'^) and thus, only local properties of the equipartition energy 
function are considered. 

Theorem 4.3. Let V he a bipartite v partition of G with rj = (f){G; V) — {u — 1) > 
and H{V] a) be the associated parameter dependent Hamiltonian. If there exists an 
equipartition {V,a), then, in the vicinity of a the set Q{V) of equipartitions forms a 
smooth 1]- dimensional submanifold ofMJ'^'^\ 

The point critical point of the equipartition energy A{V; a) on the 

manifold Q{V) if and only if it corresponds to an non- degenerate eigenfunction with 
eigenvalue A{V;a^) and nodal domains V. Furthermore, if A{V; a'^) is the n-th eigen- 
value in the spectrum, the Morse index fin of the critical point satisfies 

(35) Hn = n - Vn- 

Proof. We start by describing a parametrization of the manifold of equipartitions by 77 
independent parameters {^fc}^^^. The construction is local around an existing equipar- 
tition. Start with the partition V (see Fig. Hl^a)), and remove rj = cj) — {v — 1) edges 
which connect different nodal domains, leaving v{G) — 1 bridge edges which turn the 
partition graph G{V) into a tree (as in Fig. IHb)). The set of removed edges will be 
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(a) 



(b) 



Figure 1 . (a) A graph with a bipartite partition. Circles enclose vertices 
within the same domain, (b) The dashed edges are removed from (a) to 
generate a tree-partition. 

denoted by X{V). Every deletion of an edge e = is accompanied by the mod- 

ification of the vertex potentials at the vertices i and j by a parameter as usual. 
This partial set of 77 parameters C,e uniquely defines the equipartition energy. Indeed, 
the graph which was produced by the deletion of edges and which we denote by is 
exactly of the kind which was discussed in theorem 14.21 Since we started at an existing 
equipartition {V,a), it corresponds to an equipartition {V,P) of the graph G^. Here /3 
is a set of i^iG) — 1 entries of 5 that correspond to edges e ^ X{J^\ The 77 parameters 
^, on the other hand, are the entries of a that correspond to the edges e G X{V). 

From Theorem 14. 2 1 we conclude that the equipartition corresponds to a Courant-sharp 
eigenfunction ip^'^\G0. Under local variation of the parameters ^ the eigenfunction 

tp^'^\G^) remains Courant-sharp and induces an equipartition of the graph G^ and, 
therefore, of the graph G. Moreover, the eigenfunction is smooth as a function of the 
parameters C,- This allows us to define the parameters a (thus constructing a locally 
smooth immersion Q{V) — ?■ W^^'^^ ) in the following manner. There is one parameter for 
each edge in S{G) \ £{V). If the edge e is in the set X[V), we take ctg = ^e- Otherwise, 
we compute Ue from the eigenfunction ip^'^\G^) according to the familiar prescription, 
see equation (120|) . 

We now prove that the critical points of the equipartition energy A{V] a) correspond 
to the eigenvectors of H{G) and vice versa. The search for critical points should be 
carried out in the manifold Q{V). To do so we search for the critical points of Ai(Pi; a) 
and impose the restriction to Q{V) by introducing (z/ — 1) Lagrange multipliers {Cfc}- 
One has to search for the critical points of 



(36) 



A(P; a) = Ai(Pi; a) + G MPk, «) - Xi{Pi; a)) , 



k=l 



The sum can be written in a more concise way as 



(37) 




k=l 
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where the are hnear in the (k. For every e = only two terms in (1371) depend 

on the parameter a^, namely the terms corresponding to -P^(i) and -P^O)' '^here the 
function x{') maps a vertex to the number of the corresponding domain. Taking the 
derivative with respect to aij, we get 

dAjV; a) gAi(P^(.);«) , d\,{P^^,y,a) 
y^^^) = Cx(0 + CxO') - 



As previously, let g{Pk', a) be the first eigenvector of the domain Hamiltonian H{Pk; a), 
normalized and positive. Using first-order perturbation theory and the explicit depen- 
dence of the potential on aij, similar to what is done in theorem 13.21 the following 
equations must be satisfied at the critical point for every G S{G) \ £{V): 

(39) c^i,i)gi{P^(j,)]af - c^{'j)\gj{Px{j)\ = 0. 

ij 

We immediately conclude that the critical values of the Lagrange multipliers are 
non-negative. Form the function 

(40) f = J2±V^k9{Pk;a), 

k=l 

where the signs are to be chosen in accordance with the bipartite structure of the 
partition V. Then equation (139|) can be written as 



f-^i^ = or a, = ^. 



which describes the values of a^ at the critical points. The function / is an eigenfunc- 
tion of the Hamiltonian H{V, a'^) — in fact one belonging to the {v — l)-dimensional 
eigenspace of the degenerate lowest eigenvalue K{V\a'^). It can now be shown by ex- 
plicit calculation, see the discussion leading up to equation (fTTj) . that / is also an 
eigenfunction of the original Hamiltonian H{G). 

Conversely, starting from an eigenfunction / with the given nodal domain structure 
P, we define 

— li — \^ f2 

and check that the condition for the critical point (139|) is satisfied. 

Finally, to show (135|) . we go back to the graph and start the process of successive 
addition of edges and use (p^ at each step, since all parameters ae are negative. Thus, 
the change in the position in the spectrum An over t] additions equals the number of 
times the critical point is approached as a maximum: 

An = ^M(e). 

Since at the beginning of the process we had a Courant-sharp eigenfunction with u 
domains, its position in the spectrum was u. Thus the change An of position is from u 
to n, or, 

(41) n - i/„ = /i„, 

where is the total number of independent directions in which one approaches the 
critical point of the equipartition as a maximum; the Morse index (see [1^ Section 4.3] 
for more detail). 
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The above construction was carried out for a particular parametrization of the mani- 
fold in the space of parameters corresponding to equipartitions. However, since A(P; a) 
is analytic in the neighborhood of the critical point, the Morse index does not depend 
on the choice of coordinates. □ 
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